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We consider the disordering dynamics of an interacting binary alloy with a small admixture of 
vacancies which mediate atom-atom exchanges. Starting from a perfectly phase-segregated state, 

Qthe system is rapidly heated to a temperature in the disordered phase. A suitable disorder param- 
eter, namely, the number of broken bonds, is monitored as a function of time. Using Monte Carlo 
QO ' simulations and a coarse-grained field theory, we show that the late stages of this process exhibit 

dynamic scaling, characterized by a set of scaling functions and exponents. We discuss the uni- 
versality of these exponents and comment on some subtleties in the early stages of the disordering 
process. 
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I. INTRODUCTION 



The dynamics of binary alloys, undergoing mixing or segregation processes, is a venerable problem in the physics 
and materials science communities Q. Much interest has focused specifically on the dynamics of phase ordering 
and domain growth, following a rapid quench below the coexistence curve. Starting from an initial non-equilibrium 
configuration, the system evolves towards its final equilibrium state, giving rise to fundamental questions in both 
disciplines. Crucial to metallurgy, the domain morphology is intimately linked to structural stability while statistical 
physics focuses on the universal, self-similar aspects of its temporal evolution, manifested in characteristic growth 
exponents and dynamic scaling. 

A key ingredient in the study of these processes is the mechanism by which two particles exchange positions. From 
the viewpoint of modelling and simulation purposes, direct particle-particle exchanges obviously lead to the simplest 
codes and coarse-grained equations. Moreover, invoking universality, it is expected that the details of the microscopic 
, mechanisms do not affect large-wavelength, long-time properties, such as growth exponents and universal scaling 
t-H 1 functions. In this spirit, Kawasaki dynamics and the Cahn-Hilliard equation have been extensively used to describe 
0^ | phase ordering in binary alloys. However, in most real solids, microscopic atom-atom exchanges can be mediated by 
a variety of processes j2j , with direct exchanges playing a rather small role since steric hindrance tends to create large 
' energy barriers. The most common mechanisms involve exchanges with defects, such as vacancies or interstitial sites. 
For these reasons, alloys have often been modelled by three-state models |J whose dynamics is controlled by atom- 
vacancy exchanges, with direct exchanges being completely forbidden. Two problems in particular have attracted 
^ , considerable attention: first, the effect of vacancies on phase separation and domain growth and second, their 

role in atomic interdiffusion . A number of studies have also addressed vacancy-mediated ordering in a variety of 
antiferromagnetic alloys , as well as surface modes of unstable droplets in a stable vapor phase || . 

In this review, we will describe a third aspect of defect-mediated dynamics, namely, the "inverse" of the phase 
ordering problem. Instead of studying the growth of order in response to a sudden temperature decrease (a "quench" ), 
we focus on the disordering of a finite system, following a rapid increase in temperature. Starting with a zero- 
temperature ferromagnetic configuration, i.e., a perfectly phase-segregated system with sharp interfaces, we monitor 
how the interfaces "roughen" and how particles of one species are transported into regions dominated by the other 
species. Clearly, if the final temperature is sufficiently high, the interfaces will eventually dissolve completely, leaving 
us with a homogeneous final state. Several questions emerge quite naturally: Are there characteristic time scales on 
which the disordering takes place, and how do they depend on system size, temperature and other control parameters? 
How do local density profiles and correlation functions evolve with time? Are there any scaling regimes, and what 
are the appropriate scaling variables? How do these features respond to changes in the relative concentrations of 
vacancies and alloy components? 

In the following, we address some of these questions in a simple model for defect-mediated interface destruction 
and bulk disordering. We consider a symmetric (Ising-like) binary alloy of A and B atoms which is diluted by a very 
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small number of vacancies (defects), reflecting the minute vacancy concentrations (of the order of 10 -5 ) found in 
most real systems. Following an upquench from zero to a finite temperature T, the vacancies act as "catalysts" for 
the disordering process, exchanging with neighboring particles according to the usual energetics of the (dilute) Ising 
model. The particles themselves form a passive background whose dynamics is slaved to the defect motion. Thus, 
this system corresponds to a real material in which the characteristic time scale for vacancy diffusion is much faster 
than the ordinary bulk diffusion time. While vacancies are typically distributed uniformly in the bulk, certain defects 
may prefer to accumulate at the interfaces. Thus, the number of defects is not necessarily extensive in system size. 

While we allow for some variation in the vacancy number, we consider equal concentrations of A and B atoms. Thus, 
our work forms a natural complement to the only other study |l0| of vacancy-mediated disordering in the literature. 
There, the alloy composition is chosen highly asymmetric: 95% of A atoms versus only 5% for the B species, with 
a single vacancy. Thus, the A atoms form a matrix for a H-precipitate. The alloy is first equilibrated at a very low 
temperature, so that small clusters of B atoms are present. It is then rapidly heated to a higher temperature, and the 
number and size of B clusters are monitored. Three different scenarios are observed, depending on whether the final 
temperature is below the miscibility gap, above the miscibility gap but below T c , or above T c . In the first case, the 
precipitates remain compact. They dissolve partially at first, but then equilibrate again by coarsening. In the second 
case, the precipitates also remain compact but eventually dissolve completely, mostly through "evaporation" from 
their surfaces. In the third case, the clusters decompose rapidly ("explode") into a large number of small fragments 
which then disappear diffusively. 

The complete or partial mixing of two materials at an interface plays a key role in many physical processes, 
such as corrosion or erosion phenomena p"l| ]. We mention just two applications with huge technological potential 
for device fabrication. The first concerns nanowire etching by electron beam lithography p^ j: If a thin film of 
platinum is deposited on a silicon wafer, interdiffusion of Pt and Si produces a mixing layer. If this layer is heated 
locally by, e.g., exposure to a conventional electron beam, silicides, such as Pt2Si and PtSi, form. The unexposed 
platinum can subsequently be etched away, leaving conducting nanoscale structures behind. The second example 
concerns mesoscopic superlattice structures, consisting of alternating magnetic and nonmagnetic metallic layers. If 
adjacent magnetic layers couple antiferromagneticahV the application of a large uniform magnetic field to these 
layered structures results in giant magnetoresistance [13]. However, the performance of these devices requires precisely 
engineered layer thicknesses and interfaces, and can be significantly affected by disorder jlj], including interdiffusion 
or interfacial fluctuations. 

While motivated by these applications, our study can only form a baseline here, for further work on more realistic 
models. However, it also has some rather fundamental implications. First, it serves as a testing ground for a basic 
problem in statistical physics, namely, how a system approaches its final steady state, starting from an initial non- 
stationary configuration. A second view of our study addresses the effect of a random walker on its background 
medium. Each move of a vacancy rearranges the background atoms slightly, leaving a trail behind like a child running 
across a sandy beach. In the simplest case, the walker is purely Brownian. In our language, this corresponds to an 
upquench to infinite temperature, T — oo, where energy barriers are completely irrelevant. In this case, there is no 
feedback from the background to the local motion of the vacancy. Nevertheless, each displaced atom displays its 
own intriguing dynamics (mainly in d = 2) ]|l6|,^7). To study the collective behavior of the atoms, |L8| explored a 
lattice filled with just two species of (indistinguishable) particles. Beyond this simple case is a system heated to finite 
temperatures. Since the background affects the vacancy through the local energetics associated with the next move, 
no exact solutions are known. Instead, progress relies mainly on simulations. 

The key result of our study is the observation of three distinct temporal regimes, separated by two crossover times, 
provided the final temperature is not too close to T c . The intermediate and late stages of the disordering process 
exhibit dynamic scaling, with characteristic exponents and scaling functions that are computed analytically. For the 
system sizes considered here, a clear breakdown of these scaling forms is observed for temperatures within about 10% 
of T c . We argue that this occurs when the correlation length becomes comparable to the system size. In contrast, the 
early stage of the disordering process consists of interfacial destruction in a highly anisotropic manner jL7) . 

This article is organized as follows. We first introduce our model and define an appropriate "disorder" parameter. 
We then review the purely diffusive case, corresponding to T = oo [fTs] , Simulation results for two-dimensiononal 
systems and an exactly soluble mean-field theory in general dimension are discussed. In Section 4, we present a 
number of new results, concerning upquenches to finite temperatures where particle-particle interactions come into 
play. We conclude with some comments on the effect of external fields, upquenches to temperatures T < T c , and 
details concerning interfacial destruction in the early stages. 
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II. THE MODEL. 



In this section, we describe the dilute Ising model underlying our Monte Carlo simulations. It is defined on a 
two-dimensional (d = 2) square lattice of dimension L x L, with sites denoted by a pair of integers, r ={x,y). The 
boundary conditions are fully periodic in all directions. To model the two species of particles, black ("spin up") and 
white ("spin down", displayed as gray in the figures), and the vacancies, we introduce a spin variable a r at each site 
which can take three values: +1 (—1) if the site is occupied by a black (white) particle, and <j t = if it is empty. 
Multiple occupancy is forbidden. The numbers of black (N + ) and white (N~) particles are conserved and differ by at 
most 1. The number of vacancies, M, is much smaller: M <C N + . In fact, most of our simulations will be restricted 
to M = 1, to model the minute vacancy concentrations in real systems. In our analytic work, we will also consider a 
more general case, where M is allowed to vary with system size according to M ot L 1 . Different values of the vacancy 
number exponent 7 G [0, d] will be discussed. Clearly, the single vacancy case corresponds to 7 = 0. 

The particles and vacancies interact with one another according to a dilute Ising model: 

H[J T ] =-jJ2 W W 
(r,r'> 

with a ferromagnetic, nearest-neighbor coupling J > 0. Since the dilution is so small, the behavior of the system is 
that of the ordinary (non-dilute) two-dimensional Ising model. Thus, it has a phase transition, from a disordered to 
a phase-segregated phase, at the Onsager critical temperature, T c — 2.267... Jjks JT^]. The ground state is doubly 
degenerate. It consists of a strip of black and a strip of white particles, each filling half the system, separated by two 
planar interfaces running parallel to a lattice axis. Since the particle-particle interactions are ferromagnetic, the few 
vacancies accumulate at the interfaces. 

Next, we turn to the dynamics of the model. Only particle-hole exchanges are allowed and performed with the usual 
Metropolis p3] rate: min{l, exp(— /3A7Y)}, where A7i is the energy difference of the system before and after the jump, 
and P = l/(fcgT) is the inverse temperature. The initial configuration of the system is perfectly phase-segregated, 
with two interfaces chosen to lie along the cc-axis. The vacancies are located at random positions along one of the 
interfaces. Since this is technically a zero-temperature configuration, while the dynamics occurs at temperature T > 0, 
the vacancies will move around, disordering the interfaces and, if T > T c , dissolving them eventually. The final steady 
state is an equilibrium state of the usual two-dimensional Ising model. Therefore, many of its properties are exactly 
known @. 

In our simulations, the system sizes range from L — 30 to 60. The final temperature T, measured in units of the 
Onsager temperature T c , varied between 1.1T C and infinity. Our data are averages over 10 2 to 10 4 realizations (or 
runs) depending on the desired quality of the data. The time unit is one Monte Carlo step (MCS) which corresponds 
to M attempted particle-hole exchanges. All systems investigated equilibrate after about I0 8 MCS. 

To monitor the evolution of the system, we measure a "disorder parameter" , defined as the average number of black 
and white nearest-neighbor pairs, A(L,t), as a function of (Monte Carlo) time t. This quantity is easily related to 
the Ising energy, 

A{L,t) = ±L d + ± J {H) + 0{Li) (2) 

where (•) denotes the configurational average over runs. The correction is due to the vacancies and remains much 
smaller than the two leading terms. More detailed information is carried by the local hole and magnetization densities, 
defined respectively as 

0(r,i)= (5 CTr , ) 

^(r,t) = {a r ) (3) 

The Kronecker-<5 ensures that lattice site r is occupied by the vacancy. Non-zero values of ^(r, t) indicate an excess 
of white or black particles at lattice site r, which is obviously a sensitive measure of the disordering process. The full 
time dependence of these densities can in general only be computed within a mean-field approach. However, their 
stationary forms are easily found from exactly known properties of the two-dimensional Ising model. 

III. BROWNIAN VACANCIES: T = 00. 

In this Section, we focus on the simplest case, namely, upquenches to infinite temperature. As a consequence, the 
nearest-neighbor coupling J plays no role at all, and all attempted particle-hole exchanges are executed with unit 
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rate. Thus, the vacancies perform a Brownian random walk, regardless of their local environment. We first summarize 
our simulation results which suggest the key to the mean-field analysis, namely, a separation of time scales. We then 
introduce a suitable set of exponents and scaling forms and show that the late stages of the disordering process 
display dynamic scaling. Finally, we turn to a mean-field theory and compute these exponents and scaling functions 
analytically. 



A. Simulation Results 



To illustrate the gradual destruction of the interfaces and the disordering of the bulk, Fig. 1 shows the evolution 
of a typical configuration, in a 60 x 60 system with a single vacancy (which is represented by a white square in the 
figures). At t = 0, the vacancy is located at the interface in the center which begins to break up slowly. Eventually, 
the second interface also becomes affected, as more and more particles are transported into regions of opposite color, 
until the system finally disorders completely. For later reference, we note that the last two configurations, at 10 7 and 
10 8 MCS, are both already fully random. The disordering process is clearly reflected in the number of broken bonds: 
.4(60, t), shown in Fig. 2, increases from its minimum of O(L) for the initial configuration at t = 0, to 0(L 2 ) for the 
fully equilibrated system at t = 10 7 . One clearly distinguishes three regimes, shown schematically in the inset: an 
early regime (I), the intermediate, or scaling, regime (II), and finally a late or saturation regime (III) in which the 
system has effectively reached the steady state. Tracking the motion of the vacancy, the physical origin of these three 
regimes is easily identified. For early times (regime I), the vacancy is still localized in the vicinity of its starting point, 
far from the boundaries of the system. After a time of 0(L 2 ), however, the vacancy has explored the whole system 
and is effectively equilibrated. This marks the onset of the intermediate regime. The particle distribution is still 
strongly inhomogeneous here and does not equilibrate until the system enters the saturation regime. We emphasize 
that the second and third regimes emerge only in a finite system. In an infinite system, regime I persists for all times. 
Some aspects will be briefly summarized in the concluding remarks. 







FIG. 1. Sequence of snapshots showing the disordering process of a 60 x 60 system with T = oo. The black and gray squares 
represent the two types of particles(a = ±1) and the white square denotes the vacancy(u = 0). The configurations were 
recorded after 0, 10 4 , 10 5 , 10 6 , 10 7 and 10 8 MCS. 
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FIG. 2. Plot of the total number of broken bonds, A(L,t) vs t, for L — 60 and T = oo. It shows the emergence of an early 
regime (I), an intermediate, or scaling, regime (II), and a late or saturation regime (III). The straight reference line has slope 
0.5. 

Our key observation is that, independent of dimension, Regimes II and III exhibit dynamic scaling. To characterize 
this behavior, we introduce a set of exponents: First, the saturation value of A scales with system size according to 
limt_>oo A(L, t) = Asat(L) ~ L a . Second, in the intermediate regime, A(L,t) grows as A(L,t) ~ L a tP. Finally, the 
two crossover times ("early" and "late") scale as t e ~ L Zc , and ti ~ L z ' . Thus, the intermediate-to-late time crossover 
can be summarized by the scaling form 

A(L,t)~L a f{t/U) (4) 

with a scaling function / which satisfies fix) ~ const for x 1, and fix) ~ x 13 for x -C 1 (but large enough to fall 
within Regime II). The consistency condition A(L,ti) ~ A sa t{L) yields the scaling law 

a + ziJ3 — a. (5) 

Due to the presence of the novel exponent a, (0) differs from the familiar z{3 = a which controls surface growth in, 
e.g., the Edwards- Wilkinson Q or KPZ @ models. 

We first consider a single vacancy, in d = 2. For this case, the data give a = 2 ± 0.1, j3 = 0.5 ± 0.06, z e — 2 ± 0.2, 
zi = 4 ± 0.2, and a — ± 0.1. Clearly, these exponents satisfy the scaling law (||). Excellent data collapse is obtained 
by plotting A(L,t)/ L a versus t/L Zl , shown in Fig. 3 for a range of system sizes. The same indices are observed for 
several vacancies, provided their number remains constant as the system size L is varied, corresponding to a vacancy 
exponent 7 = 0. We note briefly that the exponent a is required here to accomplish good data collapse in dimensions 
greater than 2. 

Another natural case is 7 = d — 1, corresponding to a situation where the defects are initially "frozen" at the 
interfaces, and their number scales with the interfacial area. Such a scenario could occur in surface catalysis or in 
welding. Since our data are restricted to d — 2, we examine the case 7 = 1, i.e., M cx L. Here, the data collapse only 
if (i) curves corresponding to different system sizes are shifted by an M-dependent factor along the lnt-axis, and (ii) 
the late crossover exponent takes on the new value zi = 3 ± 0.2, while a and (3 remain unchanged. This behavior 
suggests that the exponent a depends non-trivially on 7. For 7 = 1, the data are consistent with a = 1/2, so that the 
scaling form (Q) still holds and the scaling law (||) is satisfied with the modified exponent zi. 

This concludes the discussion of our simulation results. We now turn to an analytic description of the disordering 
process. 
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FIG. 3. The scaling plot of A(L,t)/L 2 vs t/L 4 for L — 30 — 60 with T = oo. The straight reference line has slope 0.5. 

B. Mean-Field Theory 

As a first step, we summarize the exactly known results for the initial configuration and the final, fully disordered 
state. The analysis is easily performed in general dimension, for a hypercube of side length L with fully periodic 
boundary conditions. The initial value of the disorder parameter is given by the number of broken bonds across the 
two initial, flat interfaces between black and white particles: A(L,t = 0) = 2L d ~ 1 + 0{M). The saturation value, 
■Asat(L) = |L 2d / (L d — 1) follows easily from the fact that the final steady state is completely random for the Brownian 
vacancy case with T — oo. Thus, we can read off the exact result a = d. 

To proceed further, we derive a set of equations of motion, for the coarse-grained local hole and particle densities. 
Starting from the microscopic master equation, we can easily obtain evolution equations for the discrete densities, 
4>{r,t) and ^(r, t). Since there is no feedback from the particle background, the equation for <fi(r,t) is completely 
independent of ^(r, t). In contrast, the motion of the particles is slaved to that of the vacancy, so that the equation for 
ip(r 7 t) is inherently nonlinear and contains two-point functions of the form (a r 'S<j r ,o)- Here, a mean-field assumption 
is required, to truncate these averages. Finally, we take the continuum limit, by letting the lattice constant vanish at 
fixed system size L. For simplicity, we use the same notation for discrete densities and their continuous counterparts. 
The position vector r = (x\ r ...,Xd) denotes a point in the hypercube — L/2 < Xi < +L/2, i = 1,2, with volume 
V = L d . We also retain the symbol t for time, since Monte Carlo time and its coarse-grained counterpart differ only 
by a scale factor. After some suitable rescalings, we obtain the desired mean-field equations: 

8t4>(T,t) =V 2 0(r,i) 

d t ^(r, t) = 0(r, i)VV(r, t) - V(r, *)V 2 0(r, t). (6) 

The simple diffusion equation for (f> reflects the Brownian random walk of the vacancies. The magnetization density ip 
obeys a balance equation: the first term reflects a gain, provided a vacancy is initially present and a particle "diffuses" 
in from a neighboring site. The second term accounts for a loss, due to a particle jumping to a vacant nearest-neighbor 
site. Similar equations have been discussed in the context of biased diffusion of two species |25| and vacancy-mediated 
interdiffusion [[7|. Note that both equations take the form of continuity equations, due to the conservation laws on 
the particle numbers. The two densities are normalized, according to J v 4>(r, t) — M, and J v ijj(r, t) = 0. Given fully 
periodic boundary conditions, these equations have to be supplemented with appropriate initial conditions. For our 
simulations, we choose <j>(r, 0) = L h d I _ 1 5(y) where y = Xd, and ip(r, 0) = 26{y) — 1. The final state is, of course, trivial: 
0(r,oo) = (j> = M/V, and ^(r,oo) = 0. 

To solve for the hole and magnetization densities, we seek the separation of time scales, first observed in the Monte 
Carlo data, in the equations (^|). The hole density, with a diffusion coefficient of 0(1), relaxes rapidly to the final value 
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M/V . By contrast, the magnetization density is essentially slaved to <fi(r,t) and relaxes with a diffusion "coefficient" 
of 0(1/V). The hole density spreads diffusively and reaches the boundaries of the system after a time of 0(L 2 ), 
marking the end of the early regime. Thus, we identify the early crossover time t e oc L 2 and read off z e = 2. For later 
times, the vacancies are uniformly distributed over the system, so that cj> may be replaced by its stationary value, 
(f) = M- . Inserting this into the equation for tp results in a simple diffusion equation 



dt*l> = o VV 



(7) 



with a diffusion coefficient M/V. Its solution, subject to the initial and fully periodic boundary conditions, is easily 
found: 



7T ' — ' 

n=l 



sin[2ir(2n- l)y/L] 
2n-l 



-et(2n-l) 



(8) 



where e = Aw 2 (f> / L 2 . Of course, ip(r,t) depends on y only, reflecting the homogeneity of initial and boundary 
conditions in the d—1 transverse directions. We also note that, while this form is exact, it converges rapidly only for 
late times, et 3> 1, corresponding to the saturation regime. There, the profile is harmonic with a rapidly decaying 
amplitude: ip(r,t) ~ ^sin[2ny / L]e~ et . Below, we will present an equivalent form with good convergence in the 
opposite limit. 

To describe the simulation results, we need an expression for the disorder parameter A. Recalling its connection to 
the average Ising energy, Eqn (g), we invoke the continuum limit of the latter |^6|, (H) = —dJ J v r/;(r, t) 2 , which is 
correct up to surface terms of 0(1/ L d ~ v ). Thus, we obtain 



A(L,t) = - 



v 



(9) 



We note that this form is manifestly extensive. Moreover, since surface terms have been dropped, the initial value 
of A is now simply 0. 

The time evolution of this quantity follows from d§|) : 



A(L,t) = -V[l-S(2et)] 



(10) 



where 5(C) = A> J2T e- c(2 "" 1)2 /(2n - l) 2 . Since 5(0) = I and 5(oo) = 0, we verify that A does take on the 
correct initial and final values. It is also consistent with our postulated scaling form, Eqn (|^), provided we identify 
the scaling exponent a with d and the scaling variable t/ti with et. This suggests to define the late crossover time 
as ti = f/e. Collecting the dependence on system size, we identify tj cx L 2 /(f> = L d+2 -~> '. Hence, we read off 
zi = d + 2 — 7 for the scaling exponent which controls the late crossover time. Thus, et 3> 1 marks the saturation 
regime where Eqn (|J) converges well, exhibiting a simple exponential approach to the saturation value: A(L, t) = 
IV[l-£e-™ + 0(e-^)]. 

In contrast, the intermediate regime (II) corresponds to et <C 1. To capture the time dependence here, we re-express 
the infinite sum via a Poisson resummation p7| . The key advantage of this procedure is to convert a sum with rapid 
convergence in one limit (et 3> 1) into an equivalent sum with good convergence in the opposite (et -C I) limit. 
Deferring technical details to an Appendix, we introduce u m = nm/\/8et and arrive at 



A(L,t) 



2d 
n^I 2 



VV2et { I 



OO 

E 



(-iy 



(11) 



where T(», •) denotes the incomplete Gamma function. In this form, the sum over m is suppressed for small et. Thus, 
A(L,t) oc V\p2et oc L d Alt / L 2+d , yielding the remaining indices, namely (3 = h independent of dimension, and 
a = \ (d + 7 — 2) . It is straightforward to check that our exponents satisfy the scaling relation (H) . In two dimensions, 
our theory predicts zi = 4 and a = for 7 = 0, while zi — 3 and a = \ for 7 = 1, in complete agreement with the 
Monte Carlo data. 

Since the full magnetization profile contains more detailed information than the disorder parameter, one might also 
be interested in its form for early times, et <C 1. A similar Poisson resummation, of Eqn (1q), results in 



^i(r, t) — erf 



y 



2^/et 



erf 



- + erf ' 



2^/et J 



2y/et 



(12) 
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where y = 2ny/L. For et <C 1, this form still reflects the sharpness of the initial 0-function which is only gradually 
"washed out" . 

This concludes our discussion of the mean-field theory, for the Brownian vacancy case. We note briefly that our 
results remain unchanged (up to a trivial amplitude) if periodic boundary conditions are replaced by reflecting ones. 
We should also reiterate that our focus here rests on the intermediate and late regimes. A full microscopic treatment 
of the early regime is beyond the scope of this paper. For the benefit of the reader, however, a short summary is 
provided in the conclusions. Instead, we turn to the role of interactions. 

IV. THE EFFECT OF INTERACTIONS: T c < T < oo. 

When the final temperature of the system is finite, T < oo, the full Hamiltonian, Eqn ([!]), comes into play when a 
vacancy attempts to move. In particular, the defects no longer perform a simple random walk, since the jump rates now 
carry information about the local environment, i.e., the distribution of black and white particles around the originating 
and receiving site. Thus, a feedback loop between vacancy and background is established, in stark contrast to the 
case T = oo. One might expect that this has significant consequences for profiles and disorder parameter, and that 
dynamic scaling exponents might be modified. It is immediately obvious that we should distinguish final temperatures 
above criticality from those below. In the former case, the final steady-state configurations are still homogeneous, even 
if correlations become more noticeable. In the latter case, however, the equilibrium state itself is phase-separated, so 
that the only effect of the vacancies is to "soften" the initial interfaces. This aspect of vacancy-induced disordering 
should be particulary interesting in two dimensions, where Ising interfaces are known to be rough |2S| |. There, the 
vacancies must mediate both the development of an "intrinsic interfacial width" and the emergence of large-scale 
wanderings of the interface [29[] . Certainly, we expect these two phenomena to take place on drastically different time 
scales, since they are associated with local vs. global disorder. 

In this Section, we focus on final temperatures above criticality: T c < T < oo, and a single vacancy, M = 1. Starting 
with a set of typical evolution pictures, we investigate the dynamic scaling properties of the disorder parameter. Since 
the equations of motion now become highly nonlinear, we have not been able to solve them exactly. Thus, we offer 
only a few comments in conclusion. 

Fig. 4 shows the evolution of a typical configuration, in a 60 x 60 system with a single vacancy at T = 1.5T C . Note 
that the MC times are the same as in Fig. 1. Again, we observe the gradual, yet eventually complete, destruction of 
the interfaces. In comparison with Fig. 1, however, two obvious differences emerge: First, the system takes longer to 
reach the final stationary state. Second, the final configuration shows clear evidence of a finite correlation length. 
Both features are induced by the interactions. In particular, the disordering process is slowed down since the breaking 
of bonds is energetically costly. At a more quantitative level, this is documented by Fig. 5 which shows the disorder 
parameter, for a 60 x 60 system and several temperatures. One observes that the late crossover time shifts to later times 
as the temperature decreases. In contrast, the early crossover time appears to be less affected. Also, the saturation 
value of the disorder parameter decreases. However, being essentially the Ising energy, it remains extensive, with a 
T-dependent amplitude. 

It is natural to ask whether the scaling forms found in the previous section still hold. In Fig. 6, we present the 
scaling plots for the disorder parameter, i.e., we show A(L,t)/ L d plotted vs. t/L 4 . Each graph corresponds to a 
different temperature: T = 3.5T C (Fig. 6a), T = 1.5T C (Fig. 6b), and T = 1.1T C (Fig. 6c). Focusing on Fig. 6a, it is 
quite evident that the data collapse just as well for T = 3.5T C as for T = oo (cf. Fig. 3). The three temporal regimes 
are well separated and easy to distinguish. In the intermediate regime, the disorder parameter follows the t 1 / 2 power 
law, over more than three decades. Thus, the data at this temperature still obey the scaling form (Q), with the same 
set of exponents. However, for the next temperature, T — 1.5T C (Fig. 6b), the range of perfect data collapse begins 
to shrink: the curve for the smallest system size (L = 30) merges noticeably later than in Fig. 6a. Also, the power 
law no longer persists over such a wide region, and the exponent (3 becomes more difficult to determine. Still, the 
larger system sizes collapse rather well, so that there is still a well-defined scaling function even if its form appears to 
change. Finally, Fig. 6c shows the data for T = 1.1T C : here, the imposed scaling form manifestly fails to match the 
data. 
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10 6 MCS 10 7 MCS 10 8 MCS 

FIG. 4. Sequence of snapshots showing the disordering process of a 60 x 60 system with T = 1.5T C . The black and gray 
squares represent the two types of particles(a = ±1) and the white square denotes the vacancy(<r = 0). The configurations 
were recorded after 0, 10 4 , 10 5 , 10 6 , 10 7 and 10 s MCS. 

We conclude that, for temperatures not too close to the critical temperature T c , the dynamic scaling form of the 
purely diffusive case, Eqn (Q), still holds. In two dimensions, we identify the saturation exponent a = 2 and, for a 
single vacancy, the late crossover exponent z; = 4. The scaling function itself, however, develops more curvature in the 
intermediate regime, so that it becomes more difficult to extract a pure power law, t@ . One might be tempted to fit the 
data to a pure power law, resulting in a temperature- dependent exponent (3, noticeably smaller than 0.5. However, we 
believe that there is little physical support for such a drastic change in scaling properties, especially in the disordered 
phase of the Ising model. Instead, we conjecture that (3 = 1/2 is still valid but that T-dependent corrections-to-scaling 
begin to play a stronger role. In support, we checked that the data for T = 1.5T C , in the earlier part of regime II, can be 
well fitted by a power law with an exponential correction, i.e., A(L,t)/L 2 cx t 1 / 2 {l + Ci(T) exp[— C2(T)i/L 4 ]}. This 
form is gleaned from Eqn by keeping the leading and first sub- leading term but allowing for two temperature- 
dependent fit parameters C\(T) and C2(T). Of course, we should anticipate significant changes in the scaling behavior 
upon entering the critical region. Since 1.1T C is "at the doorstep" of the latter, some precursors such as an increased 
correlation length begin to make their influence felt. 

In the following, we comment briefly on some analytic results. First, we can easily check that the saturation values, 
■Asat{L), are consistent with exact results for the two-dimensional Ising value. Table 1 shows this comparison for a 
60 x 60 system with a single vacancy, at several temperatures. Ag a f is the measured saturation value of the disorder 
parameter. A T s ^ t is calculated on the basis of Eqn (0), where (Tt) is the exact bulk value for the average energy of 
the usual Ising model pi] ]. The effect of the vacancy has been neglected. The agreement is very good, within the 
statistical errors of our data. The discrepancies are of course largest for the last column, with T closest to T c . 
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icP 10 2 10 4 10 6 10 8 
t (MCS) 

FIG. 5. A(L,t) vs t for L = 60 with T = oo, 3.5T C , 1.5T C and 1.1T C . 

Finally, since (H), and therefore *4 sat (L), are extensive, we still have an exact value for the saturation exponent 
a = d. 

The time-dependence of the disorder parameter poses a much more difficult problem. It is of course possible to 
obtain the temperature-dependent generalization of Eqns (||) , by averaging the microscopic master equation, followed 
by a mean-field approximation and a naive continuum limit. For a dynamics slightly different from ours, this procedure 
was followed in j^]. However, these equations are no longer easily soluble: First, both of them are inherently nonlinear, 
due to interactions. Second, the equation for the local hole density, 0(r, t), no longer decouples from the magnetization 
ip(r,t), since the background feeds back into the motion of the vacancy. Thus, it is not apparent how a clean time 
scale separation could emerge from these equations. As a consequence, we are reduced to numerical solutions or 
approximate techniques, such as a high temperature expansion. Preliminary results, based on the latter, provide 
some support for our conjecture that (3 = 1/2 is still valid and that larger corrections-to-scaling are responsible for 
the curvature in the scaling function pofl . 

To summarize this section briefly, our data indicate that, for temperatures not too close to the critical temperature 
T c , the dynamic scaling form of the purely diffusive case still holds, even though corrections-to-scaling become more 
noticeable. In contrast, for T near T c we observe a clear breakdown of these scaling forms. 



TABLE I. Comparison of exact and measured values of the saturated disorder parameter for several temperatures. 



T/T c 


OO 


3.5 


1.5 


1.1 




3600.0 


3134.6 


2367.0 


1602.9 




3598.6 


3134.0 


2368.9 


1615.7 



10 




10" 8 10" 6 10" 4 10" 2 icP 10 2 
c) t/L 4 

FIG. 6. The scaling plot of A(L,t)/L 2 vs t/L 4 for L = 30 - 60 with (a) T = 3.5T C , (b) T = 1.5T C and (c) T = 1.1T C . The 
straight reference lines have slope 0.5. 

V. CONCLUDING REMARKS 

To summarize, we have analyzed the vacancy-driven disordering process of an initially phase segregated binary 
system, which is rapidly heated to a temperature T above criticality. The number of vacancies scales with system 
size as L 1 , where the vacancy exponent 7 is a parameter. To quantify the evolution of the system, we measure 
the number of broken bonds, A(L,t), as a function of time. This "disorder parameter" allows us to identify three 
temporal regimes, distinguished by the distribution of the vacancies through the system. Our key result is that the 
late stages of this process exhibit dynamic scaling. A set of exponents {z e , zi, a, a; /?} can be defined, characterizing, 
respectively, the system size dependence of two crossover times, the final saturation value of the disorder parameter, 
and its amplitude in the intermediate regime. The temporal growth of A during the latter regime is captured by the 
exponent /3. If the final temperature is infinite, i.e., T = 00, the motion of the vacancies is a simple Brownian random 
walk, and all indices can be computed analytically, in excellent agreement with the data. In particular, we find that 
the typical time scale U ~ L Zl , which controls the crossover between increasing disorder and saturation, is set by 
zi = 2 + d — 7, and thus depends explicitly on both the space and vacancy dimensionalities, d and 7. Measurements of 
Zi can therefore provide information about the vacancy distribution in a sample. In the most familiar case, standard 
vacancy diffusion in solids, the number of vacancies is extensive (7 = d), so that the well-known result zi — 2 is 
reproduced In contrast, in two dimensions we find zi = 4 for a single vacancy (7 = 0) and zi = 3 if the defects 
are generated at the initial interface. For finite temperature T > 1.5T C , we observe that these scaling forms still 
hold. Even though interparticle interactions now play a role, correlations in the system are still short-ranged, so 
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that the vacancy still performs a random walk if viewed on a length scale which exceeds the correlation length, £(T). 
As T approaches T c , however, £(T) reaches O(L) so that the simple random walk scenario must break down. In 
our simulations, significant deviations from the high-temperature scaling form appear already for T = 1.1T C , where 
£(T) ~ 6, in units of the lattice spacing 

All the investigations above are focused on upquenches to the disordered state. A natural extension of this project 
is to set the final temperature below criticality. Needless to say, new questions and interesting phenomena can be 
expected. In particular, the final equilibrium state is still phase separated, so that all the interfacial properties come 
into play. With a single vacancy in a finite system, there may be a further separation of times scales. We expect that, 
sooner after the early stage, the "intrinsic" density profile of the interface can be built. At this point, the interfacial 
width, w, is most likely controlled by £. On the other hand, at sufficiently long times, the capillary waves will surely 
make their presence felt. In the d — 2 case, the interface is always rough, so that we can always expect the latter 
crossover to occur. How, and if, w(^,L;t) scales will be of great interest. Besides dynamics and timescales, we could 
study the equilibrium probability profile of the vacancy, i.e., where the vacancy spends most of its time. In case more 
than one vacancy is present, this profile should map into a density profile for the vacancies. In either case, we expect 
the vacancies to be trapped, to a greater or lesser extent, at the interface. Now, if more and more vacancies are 
added to system, the interface profile should be altered, since a preponderance of vacancies may significantly modify 
the interfacial energy. For d — 3 Ising models, there are further interesting phenomena, associated with roughening 
transitions. Given that interfacial energies should depend on vacancy concentrations, how the locations (if not the 
nature) of such transitions are affected is a natural question. Finally, it would be extremely interesting to study the 
dynamic content of these systems, such as scaling properties when the upquench is set at a roughening temperature. 

We conclude with a few remarks on related problems. It is obviously interesting to test the range of universality of 
the high-temperature scaling forms, with respect to changes in the boundary conditions and the microscopic dynamics. 
Returning to T = oo, we have investigated several cases. First, replacing periodic by reflecting boundary conditions 
changes only trivial amplitudes. A more significant modification is the addition of a bias E, such as a gravitational 
or electric field pl| : if all particles carry equal charge while the vacancy is neutral, vacancy moves along a selected 
direction will be suppressed. In this case, the boundary conditions determine the nature of the resulting steady state. 
For reflecting boundary conditions, it is described by a simple Boltzmann factor, incorporating a gravitational or an 
electrostatic potential. Here, a new scaling variable, A = LE, emerges which characterizes the potential difference 
between opposing boundaries. In contrast, in the periodic case the steady state is a non- equilibrium one. Moreover, 
it is translation- invariant, so that all dependence on E disappears in a suitably chosen co-moving frame! Remarkably, 
however, these boundary-induced differences appear only in the scaling functions, leaving us with a universal set of 
scaling exponents {z e , zi, a, a; (3} [|TJ. 

Finally, we should add a few remarks on the early regime. As pointed out above, for an infinite system with a 
single vacancy, this is the only regime. First, notice that only the phenomenon of interfacial destruction is relevant, 
since bulk disordering would take infinitely long times. Next, considering the effective range of a Brownian walker, 
it is clear that the size of the disordered region grows at most with t 1 ! 2 . In more detail, we see that there are two 
distinct types of disorder. Loosely, we will label them "transverse" and "longitudinal." The first type is associated 
with a monolayer-like disorder from the initially sharp interface, spreading from where the vacancy began its journey. 
In particular, each time the vacancy crosses the plane at a new location, a particle is moved into the "wrong" phase. 
Thus, we may expect transverse disorder to spread like t 1 ^ 2 . The second type is a measure of the width of the 
interfacial profile. For this to occur, the vacancy must "carry" e.g., a black particle deeper and deeper into the sea 
of white ones. Thus, the return probability of the walker plays a crucial role. Only in d = 2 is this probability both 
relevant and non-trivial. Applying the well-known results of random walks |2^ | , we find that "longitudinal" disorder 
grows as (lni) 1 ^ 2 |l7j]. 

Even though our model is very simple, it forms the basis for the description of a large variety of related systems. 
Moreover, it is truly gratifying that considerable analytic progress is possible for a problem that is both nonlinear 
and time- dependent. Work is in progress to probe the effect of fluctuations on the structural decomposition process 
at a more microscopic level p2[ . 
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APPENDIX: POISSON RESUMMATION 



In this Appendix, we add a few details of the Poisson resummation procedure. Or starting point is Eqn jicj). 
Dropping the constant prefactor, we focus on the expression in the brackets, with £ = et: 



f(()El-^^M 2 /(2l! _ lf 



(13) 



This sum converges rapidly in the limit £ — > oo. We seek a form which converges equally well in the limit £ — > 0. We 
first recast this sum in terms of a much simpler one, namely, 



"(An) 2 



(14) 



n=l 



We begin by noting that Y^T l/(2n — l) 2 = 7r 2 /8 so that 



7T 2 ^ (2n-l) 

n=l v ' 

a /- 1 x 



-C(2n-l) 2 s 



(15) 



It is easy to recognize the sum in this expression as the "odd" terms in S(X), with A = \/s(. This suggests that we 
should recast 5(A) in the form 



5(A) = J2 ?~ (2Xnf + E e- (A(2 ™" 1))2 

ri— 1 n—1 

= 5 e (A) + 5 (A) 

where the "even" terms just reproduce 5(A), via 5 e (A) = 5(2A). Summarizing so far, we obtain 



F(C) = ^C ! da S (^C~s 
n Jo 



= -?C/ ds 



5(VCs)-5(2^Cs) 



Thus, it is sufficient for us to resum 5(A). We recall the Poisson resummation formula 
defined on < x < oo with lim^^oo f(x) = and Jl dx f(x) finite: 



1 oo 



(16) 

for a function / (x) , 
(17) 



where / is the Fourier cosine transform of /: 

/•oo 

/(a) = 2 / g?x /(x) cos(ax) . 
Jo 

From Eqns ([li]), we can immediately read off /(x) — exp(— x 2 ) for our case whence 

/(a) = V^exp(-a 2 /4) . 

Thus, we obtain 



5(A) = -i + A- 1 (^ + ^E' 



.-(m7r/A) 2 



(18) 
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Clearly, the original form of S(X), Eqn ( p"i| ) converges rapidly for large A, while the alternate form presented here 
converges rapidly for A — > 0. Inserting our result into Eqn ( |l6|) and performing the integration results in 



no 



r 3/2 



V^h + 2g(-l)» 



(19) 



where u m = 7rm/v8ei. This leads immediately to our result for the disorder parameter, Eqn (|l] 

The Poisson resummation of the magnetization profile, Eqn (j^), is a little more involved, but follows an analogous 
series of steps. In this case, the problem can be reduced to resumming 



Z J in 



71=1 



where y = lixyjh, since we can rewrite 



TT 



E(y,C)--S(2y,40 



The resummation leads to 



E(y,0 = -| + |erf 



2\/C 



IE 



erf 



y + 27rm 
2VC 



y — 2irm 
2V( 



so that we finally arrive at Eqn ( |l2|) 
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